Magnetic correlations on the full chains of Ortho-II YBa2Cu306.5 
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We propose that the NMR Une shape on the chain Cu in the stoichiometric high-Tc supercon- 
ductor Ortho-II YBa2Cu306.5 is determined by the magnetization induced on Cu near O vacancies, 
due to strong magnetic correlations in the chains. An unrestricted Hartree-Fock calculation of a 
coupled chain-plane Hubbard model with nearest-neighbor d-wave pairing interaction shows that 
the broadening of NMR lines is consistent with disorder-induced magnetization at low temperatures. 
In addition, we give a possible explanation of the anomalous bimodal line shape observed at high 
temperatures in terms of nonuniform Cu valence in the chains. The proximity between chains and 
CuO plane induces anisotropic magnetization on the planar Cu, and broadens the plane NMR lines 
in accordance with that of the chain lines, in agreement with experiment. We discuss implications 
of the model for other experiments on underdoped YBCO. 

(November 27, 2008) 
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I. INTRODUCTION 

In materials with strong electronic correlations like the 
cuprates, the response of the system to disorder can be 
quite different from weakly interacting metals. One well- 
known anomalous effect observed in disordered cuprates 
is the extended staggered magnetic droplets carrying a 
net moment which form around nonmagnetic impurities 
introduced on the Cu02 plane. This picture was pro- 
posed on the basis of extensive NMR studies, and can be 
understood within various theories of a potential scat- 
terer in a correlated host metali. The local susceptibility 
of these magnetic droplets manifests a Curie- Weiss be- 
havior with a Weiss temperature which drops rapidly as 
the materials are underdoped^. While this picture was 
established for Zn, Li, and strong scatterers in the Cu02 
plane created by electron irradiation, the magnetic re- 
sponse to the weaker potentials created by out-of-plane 
dopants is also expected to be enhanced, and has been 
claimed to be responsible for glassy behavior observed at 
low T and doping in intrinsically disordered systems such 
as Lai_a;Sr2:Cu04and Bi2Sr2CaCu208^. In such materi- 
als, it is generally believed that quasiparticles moving in 
the Cu02 plane are subjected to weak, extended poten- 
tials caused by the dopants, giving rise to small momen- 
tum transfer scattering in the electronic transport^i^i^. 

In the YBa2Cu307_5(YBCO) family where most of the 
impurity induced magnetization effects were studied, the 
system is doped in a different manner, which depends on 
the distribution of O atoms in the chain layers. Care- 
ful annealing under uniaxial pressure has been shown 
to give high quality single crystals where the chain O 
is ordered with long correlation lengths^. The so-called 
Ortho-I, Ortho-II, etc. YBCO crystals formed by these 
methods are the only known doped, ordered, stoichiomet- 
ric high-Tc materials, ideal for studying the underlying 
intrinsic physics of the cuprates without complications 
from dopant disorder. It is then interesting to explore 



the effects of the few isolated defects which remain in 
these crystals, expected to be O vacancies in otherwise 
full CuO chains^. 

The one dimensional nature of the chains them- 
selves induces, in the filled chain compound, a — b 
anisotropy which has been observed in dc-^'-- and opti- 
cal conductivityj^- penetration depth measurements ; 
and is expected to affect the vortex core structurei^. 
Several authors have proposed that the anisotropy in 
penetration depth can be explained by assuming a 
metallic chain that couples to the plane via interlayer 
hoppingiiii^iiSii^SiiSi. The superconductivity observed in 
the chains is assumed to be due to the proximity coupling 
to the plane, where pairing occurs. While initially quan- 
titative details of penetration depth disagreed with the 
proximity models, Atkinson has argued that accounting 
for disorder can remove these discrepancies^^. 

The presence of disorder has been shown to modify 
the local electronic structure of the chain layer. Scan- 
ning tunnelling spectroscopy (STS) measurements on the 
chain layer of optimally doped YBCO show a clear con- 
ductance modulatiou j^^i^^i^^ whose wavelength displays a 
strong energy dependences^! Such a dispersion in a quan- 
tity directly related to the local density of states (LDOS) 
suggests a Friedel-type oscillation description of this phe- 
nomenon, in contrast to an explanation in terms of a 
charge density wave(CDW), which usually has a fixed 
wave vector. In addition, resonance peaks that appear 
in pairs are found at small frequencies, reminiscent of 
resonant magnetic impurity states in an unconventional 
superconductor. Morr and Balatsky^^ argued that such 
features could be caused by magnetic impurities in the 
chain layer, without discussing the origin of such defects. 

In this paper we study the YBCO compound with a 
particular out-of-plane oxygen distribution, namely the 
Ortho-II structure of YBC06.5 where every other chain 
is filled, and show that the existence of magnetic cor- 
relations and their interplay with defects on the chain 
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can explain a series of NMR experiments by Yamani 
et a&^. We show in Section II that the magnetiza- 
tion resulting from the Friedel oscillations of an uncorre- 
lated system near the chain ends gives NMR linewidths 
much smaller than what has been reported. The line 
width in the uncorrelated case also remains temperature- 
independent, contradicting the experimentally observed 
broadening with decreasing temperature. Including mag- 
netic correlations in the chain layer within weak coupling 
mean field theory, we show that both the significantly 
larger linewidth and temperature dependence can be ac- 
counted for. This demonstrates the importance of mag- 
netic correlations in the chains of YBCO systems, which 
is frequently neglected. 

One of the unusual features of the lineshape observed 
on the chain Cu's is a low-frequency satellite peak. We 
show here in Section III that satellite features do not arise 
naturally due to O chain vacancies, but must be related 
to a set of sites in the material which have zero spin shift. 
One possibility we discuss in some detail is that some of 
the chain Cu's are in the Cu'^-I- configuration. 

In Section IV, we consider the coupling of the chain to 
the Cu02 plane, and show that the major effect of the 
chain-plane coupling is to induce an anisotropic magne- 
tization pattern on the plane, where the chains imprints 
their ID correlation length onto the plane. The broaden- 
ing of plane NMR lines in accordance with that of chain 
NMR lines is reproduced in the proposed model, con- 
sistent with what is observed in the correlation between 
chain and plane linewidth in experiment. 



II. MODEL HAMILTONIAN AND NMR 
SPECTRUM 

The proper choice of a model for the Ortho-II sys- 
tem relies on its unique lattice structure. Due to miss- 
ing oxygens on the empty chain, the Cu( IE) (chain Cu 
in "empty" chain) sites are highly localized and do not 
directly couple to the full chain; we therefore assume 
they do not affect the magnetic properties on the full 
chain and drop these degrees of freedom. The resulting 
minimum effective model contains a single square lat- 
tice of Cu(2E) (planar Cu above empty chain site) and 
Cu(2F) (planar Cu above full chain site) where supercon- 
ductivity occurs, couples to evenly spaced one dimen- 
sional chains via interlayer hopping, as shown in Fig. [21 
The full Hamiltonian consists of 



H — Hp + He + Hinter + Hijyip , (1) 

where Hp describes the planar hopping, pairing and mag- 
netic correlations 



= E + E(^s - A^')"- (2) 

ijcr i<7 

+ T.U'nf.nfi-T.^n^^f^'i, (3) 
and on the chain layer 

He = Y.-t%dlJ^+sa + h.c.) (4) 

xa 

+ ^(e^-M^R. + E^'^^T^^i ' (5) 

xa X 

with the coupling between them 

Hinter = E "^"^ (^L^xa + h.C.) . (6) 
{ix)(7 

We denote by Cic{dxa) the plane (chain) operator located 
at site i(x), and tfjif^) the planar(chain) hopping, with on 
site Coulomb repulsion UP{U'^) on the plane(chain). The 
chain hopping is only between nearest neighbors, while 
planar hopping t^^ contains both nearest and next 
nearest neighbor sites tj^^^ which is necessary to produce 
the Fermi surface of the YBCO system. The short range 
attractive interaction V accounts for the d-wave pairing 
in the plane, and the effect of magnetic field is included 
in the Zeeman term e^^"^ = gfiBB/2. The proper choice 
of all parameters except the U's relies on three criteria: 
(l)the Fermi surface of ortho-II calculated from density 
functional theory(DFT)^i^ is reproduced in this sim- 
ple 3-site model; (2)the homogeneous magnetization on 
the full chains is consistent with the value indicated by 
the main line of Cu(lF) NMR spectrum, and (3) a crit- 
ical temperature of Tc = 60K is obtained. Regarding 
the first point. Fig. [3 shows the Fermi surface produced 
by the present model in the normal state with V = 0. 
Although one does not expect the splitting of the plane 
band since the chains are connected to a single plane in 
this minimum model, the topology and the folding of the 
Brillouin zone is well reproduced. Notice that the chain 
band is almost parallel to the zone boundary with a very 
small intercept ky ~ 7r/4, indicating the small filling of 
the conduction electrons on the chain, which is a natural 
result of fitting the complex Fermi surface by this simple 
tight binding model. We therefore choose the chemical 
potential such that the average chain filling is 0.25, and 
the average plane filling is 0.9. Since the main line of 
the NMR spectrum is given by the sites far from the de- 
fects if one assumes a dilute impurity concentration, we 
found that = 2 gives a Knight shift and magnetization 
consistent with the value corresponding to the measured 
Cu(lF) main line, which then fixes the absolute scale of 
the band width for the chain. Finally, Tc = 60K in the 
presence of the chosen interlayer hopping fixes the pair- 
ing interaction at = 0.5, and the absolute scale of our 
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energy unit fP^^ is chosen such that the magnetic field 
scale and temperature dependence of the NMR lines are 
consistent with experiments, which gives iJin = ISOmeV. 




FIG. 1: Schematic of chain-plane proximity model for Ortho- 

II YBCO FIG. 2: (color online). Fermi surface produced by model de- 

scribed in Fig. [T]with parameters(f^„ = l,t^„n = — 0.4, t'^ = 
2,r = 0.1, fiP = -1.1, = -3.25), and we take f^n = 
ISOmeV. 



A Hartree-Fock-Gorkovmean field factorization is then 
applied to the above Hamiltonian, with order parameters 
defined as 



self-consistent equations 
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A Bogoliubov transformation is applied for both chain 
and plane operators 
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(8) 



kV^t) = El<.TlV(i?nT) + K,Tl'(l-/(^"i)). 

n>0 

(4ica) = El<al'/(^"i) + l<al'(i-/(^"T))> 

n>0 
n>0 

= El<-il'/(^»i) + K,.il'(l-/(£^nT)) , 
n>0 



n>0 



n>0 

+ <*M<,x+5lfiEnl) 



(9) 



where n > indicates only eigenstates with positive 
eigenenergies are included in the summation. Note that 
by estimating the concentration of terminal Cu, the aver- 
age chain length that produces reported NMR spectrum 
was reported to be around 1206, where b is the chain 
lattice constant by Yamani et al. We therefore simulate 
such a system by a rectangular plane with size 12 x 121, 
where each odd chain allows for hopping along its length 
of 121 sites directed along the x-direction. 



and the order parameters are determined by the following 



The resonance frequency Ui at site i is calculated by 
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converting magnetization to Knight shift Kspv, 
Sx 



^ 9^^B 



Xx = 



A 



hf 



(10) 



which yields, when combined with the orbital shift Korb, 
the local resonance ft'equency 



(11) 



For Cu(lF) NMR lines, we use Ahf = SOkOe, g = 2, 
Korb = 1-2%, and -f/2n = 11.285MHz/T. The histogram 
of collecting Vi on all Cu(lF) sites is then artificially 
broadened by a Lorentzian to give a continuous spectrum 



(12) 



where N{vi) is number of sites that have frequency Vi^ 
with 77 = 0.04 and R is the proper numerical factor that 
normalizes the area under the curve. 
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FIG. 3: (color online) Magnetization (left) and NMR spec- 
trum(right) produced by reduction of a single hopping in 
an isolated ID chain (located at the middle of the chain as 
shown), equivalent to Hamiltonian H chain + H^^^, in an ex- 
ternal field 7.7T at 50K. The enhancement of magnetization 
is consistent with the broadening of NMR linewidth as is 
increased. 

We now examine different impurity models that are 
relevant to the NMR experiments. The most straight- 
forward proposal for Himp comes from the observation 
that the chain hopping is a result of orbital overlap be- 
tween Cu(lF) and its two adjacent oxygens, such that 
abrupt termination of the full chain due to missing oxy- 
gens should dramatically reduce the hopping at the chain 
ends. A single isolated reduction of is modelled by an 
additional term in the Hamiltonian 



(13) 



with a further assumption of complete elimination Sf^ = 
—f^ and one impurity per chain. The Cu(lF) line caused 
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FIG. 4: (color online) Temperature dependence of an NMR 

lines for an isolated chain with sparse bond impurities, 

H chain + ^^j^mp ' uucorrelatcd chain(left), and the cor- 
related case(right). 



by a single hopping reduction in the chain layer is shown 
in Fig. [3l where we highlight the effect of magnetic corre- 
lation with the model containing only H chain + H^^^, 
to isolate the chain from proximity to the plane in the 
first analysis. 

The parameter range explored for is such that it is 
close but smaller than the critical value beyond which the 
system flows into long range magnetic order in the pres- 
ence of impurities, and the induced moment remains in 
the paramagnetic region. Fig. [3] shows that the broaden- 
ing of the Cu(lF) line coincides with the enhancement of 
real space magnetization as is increased, while the 
linewidth for the uncorrelated case C/"^ = is signifi- 
cantly smaller. The periodicity of magnetic oscillation, 
as well as the density modulation for electrons of both 
spin species, is of the order of 10 lattice constants from 
the defect, consistent with the small Fermi momentum of 
the chain band ky i:/4, resulting from fitting the Fermi 
surface with the 3-site model. Moreover, the tempera- 
ture dependence of the linewidth shown in Fig. 3] indi- 
cates that only the correlated case can properly account 
for the experimentally observed broadening as tempera- 
ture is lowered, while the width in the uncorrelated case 
remains constant. The linewidth plus its temperature 
dependence therefore proves the existence of magnetic 
correlations on the full chains of the ortho-II YBC06.5. 
Our best fit to the line width gives — 1.65, and will be 
the value used for the rest of the paper unless otherwise 
specified. We make detailed comparison with experiment 
in Section IV. 



III. SATELLITE FEATURE IN CHAIN Cu NMR 
LINE 

Although the Cu(lF) linewidth and temperature de- 
pendence is well described by impurity induced magne- 
tization due to ^ one important feature regarding the 
Cu(lF) line shape seems to be outside of this scheme. 
At temperature higher than TOAT or so, a satellite peak 
with significant weight gradually develops at frequency 
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slightly lower than the main line^" . Our calculations for a 
single defect or chain end indicate that magnetization de- 
cays smoothly away from the chain ends, however, lead- 
ing to no particular satellite feature in the histogram of a 
single defect. We therefore propose an alternate explana- 
tion for the satellite peak observed. Two aspects of this 
high temperature satellite appear to us to be important: 

(1) The weight of the peak is about 10% of the whole spec- 
trum, and (2)the position of the peak remains roughly 
the same at all temperatures. Aspect (1) then implies 
that roughly 10% of the sample has the same magneti- 
zation value, which results in this satellite peak. Aspect 

(2) suggests that this magnetization remains constant at 
high temperature. Since increasing temperature should 
continuously reduce any finite magnetization developed 
due to the correlation effect, it is reasonable to assume 
this observed constant magnetization is zero. In fact, the 
frequency of this satellite peak corresponds very closely 
to an NMR shift v with zero Knight shift Kspin = 0. 
We therefore propose that the high temperature satellite 
peak is due to a section of adjacent chain Cu which have 
their conduction electrons missing, i.e. Cu"*" ions. The 
valence changes may be correlated with the chain ends 
due to the change of local chemical environment which 
prohibits electrons from populating these sites, but need 
not be. Indeed the valence of the Cu(l) ions has been 
controversial, and our analysis suggests that the valence 
may be distributed quite inhomogeneously. Eliminating 
these adjacent Cu(lF) Cu in the one band model gives 
the following perturbation 



H. 



(2) 



1</<L,(T 

Y ~5f{cldi, + h.c.) 

l<l<L,i,<7 

1<1<L,<7 



(14) 



where Sf^ = —f^ eliminates the hopping between these 
chain sites, and Sf^ = — f eliminates the interlayer hop- 
ping that connects the plane sites with these chain sites. 
A strong impurity potential Uimp = 100 is introduced 
to artificially project out electrons on these sites. Notice 
that the summation over impurity position in Eq. (jl4p is 
arbitrarily restricted between a section of adjacent sites 
l<l <L with L = 15. 

The effect of removing these Cu spins is first examined 
in a single isolated chain, where the elimination of 5V 
in Eq. is omitted since all couplings to the plane 

are ignored in this approximation. In comparison with 



H, 



(1) 



randomly distributed missing oxygens Hchai,, . --„np' 

Cu(lF) lines given by considering Hcham + HlJ^ show 
a clear asymmetric line shape with significantly more 
weight at lower frequency, as shown in Fig. [5] The 
zero magnetization peak appears at all temperatures, and 
is specially noticeable at high temperature as the main 
line narrows. Broadening of the main line smears out 



the satellite peak, which may explain its apparent dis- 
appearnce at low temperature in the work of Yamani et 
al.-''. A typical real space magnetization pattern of the 
eliminated Cu(lF) model is shown in Fig. [SI where one 
can clearly identify each sector of the magnetization with 
the corresponding features in the NMR spectrum. 
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FIG. 5: (color online) Comparison of experimental ®''Cu(lF) 
NMR line of Ortho-II YBC06.5 at external field 7.7T along a 
direction(upper left), with three models studied: isolated ID 
chain with random impurities Hchain + 7^1^^ (upper right), 
isolated ID chain with a section of consecutive Cu erased 
Hchain + H\^^{\awex left), chain-plane coupled system with a 

(2) 



section of chain Cu erased Hp 
(lower right). 



+ H, 



IV. EFFECT OF CHAIN-PLANE COUPLING 

The conclusion from examining the single isolated 
chain does not alter as the chains are coupled to the 
plane; in other words, we still see (l)broadening of spec- 
trum at low temperature, (2)asymmetric spectral weight, 
and (3)the high temperature satellite peak, as shown in 
Fig. [51 The interlayer coupling influences the pairing 
and magnetic state of the plane, as we show below, but 
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biguously proves that, when systems with different di- 
mensions are coupled, the response of lower dimension 
can drive that of the higher dimension, consistent with 
our expectation that correlation effect becomes more sig- 
nificant as dimensionality is lowered. 
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FIG. 6: (color online). Real space magnetization on one of the 

chains in the Cu(lF) model Hpiane + Hchain + Hi„ter + H^^^, 
at temperature 7QK. The spinless sites, here assumed to be 
near the chain end, give zero magnetization and correspond 
to the satellite peak on the NMR spectrum. The small mag- 
netization region away from the chain end corresponds to the 
main line, while the large oscillation sector close to the spin- 
less sites broadens the main line. 



as far as the chain magnetization is concerned, the inter- 
layer coupling does not qualitatively affect its magnitude 
or spatial distribution. 

Our basic finding here is that the chain magnetization, 
which has a strongly ID character, imprints itself on the 
magnetic response of the plane to which it is coupled, 
even though the coupling itself is only via interplanar 
hopping (there are no explicit chain-plane exchange in- 
teractions included in the Hamiltonian). Disorder in the 
plane will of course produce a response with local 2D 
(4fold) symmetry, as in e.g. the occasional Cu vacancy 
which will produce a unitary scatterer. We neglect these 
effects in the current model, and focus only on the chain 
disorder which is expected to dominate in the ultraclean 
Ortho-II crystals. 

To demonstrate this pronounced planar ID anisotropy, 
as well as to compare with the planar unitary scatterers, 
we employed again the reduction of a single hopping on 
the chain iff^p as discussed in Section II, but this time 
consider the full chain-plane system Hpiane + Hchain + 

Hinter + h\^^. Fig. [Tj shows the resulting magnetiza- FIG. 7: (color online) Magnetization induced (a)on the chain 
tion pattern induced on both chain layer and the plane, 1.5,^0^ ^^a +v,t, k,. o r,^T,,^;„„ ^^^0+ h-(i) 

where one again sees the long wave length oscillation as- 
sociated with the small Fermi momentum of the chain 
band, similar to the result of Fig. |3| In addition, the 
proximity between the two layers inherits this small chain 
Fermi momentum to the plane, results in a small planar 
induced magnetization that has wave length similar to 
that of the chain. The planar magnetization is obviously 
directed along the chain direction, which confirms our 
suspicion that out-of-plane oxygen disorder in YBCO can 
induce anisotropic density and magnetic modulations on 
the plane. On the other hand, numerics show that if one 
places a unitary scatter on the plane, then the anisotropy 
of its induced magnetization is much smaller even when 
interlayer coupling is present. Such a hierarchy unam- 
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layer and (b)on the plane by a hopping reduction defect Hi 
located on a single chain, together with their corresponding 
false color plots (c) and (d) which manifest the similarity be- 
tween their periodicity. Size of the plane is 72 x 24, at tem- 
perature 7QK. 



Finally, we address the issue of correlation between 
broadening of Cu(lF) and Cu(2E/F) lines. Experi- 
mentally the identification of each line on the complete 
NMR spectrum is made via comparison with YBC06 
and YBC07 spectra, and the Knight shift can be ex- 
tracted by subtracting the orbital contribution Korb as- 
sociated with each copper species^^S The reverse process 
provides the recipe of recovering the NMR spectrum in 
the present model, with each line calculated by coUec- 
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tion of magnetization and Eq. ((TT|) . The value of Korb 
at external field B//h is taken from Takigawa et alM. 
In Fig. [HI we show that the broadening of extracted 
Cu(2F) and Cu(2E) lines as temperature is lowered fol- 
lows that of Cu(lF), with the slope of Cu(2E) smaller 
than Cu(2F). Increasing the interlayer hopping increases 
the Cu(lF) and Cu(lE) linewidth, indicating that the 
magnetization on the plane comes from its proximity to 
the chains. Within the range of interlayer hopping that 
retains the Fermi surface topology, the linewidth correla- 
tion between the planar Cu and the chain Cu is roughly 
linear, similar to the conclusion found experimentally.'^*' 
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FIG. 8: (color online) The standard deviation of Cu(lF) 
line versus the Lorentzian linewidth of the (a)Cu(2F) 
and (b)Cu(2E) line, extracted from the chain-plane model 
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Note that due to the 



plane i ^-^ cnairL i -'-^■i/tLtv i ^-^imp' 

high temperature satellite peak and the complex line shape 
of Cu(lF), as shown in Fig. [S] we use the standard devia- 
tion of the Cu(lF) line to represent its linewidth. Data taken 
at system size 101 x 10 with eliminated sites L = 11 in im- 
purity model Eq. (|14p . and the temperature of each point 
corresponds to the same scale in Fig. [5] 



V. CONCLUSIONS 

In summary, we have shown that the broadening of 
Cu(lF) NMR hues of ortho-II YBCO in the supercon- 
ducting state at low temperatures cannot be explained 
unless one assumes the existence of intermediate strength 
magnetic correlations on the chains. We introduced a 
model with correlations described by Hubbard interac- 
tions U on both the plane and chains, together with a d- 



wave pairing interaction, both of which were treated in an 
unrestricted mean field approximation. In the presence 
of an impurity, here assumed to be an oxygen vacancy in 
the nearly full chain, we find that a staggered incommen- 
surate magnetization oscillation grows as temperature is 
lowered (see Ref. [l]); many such defects give rise to a 
broadening of NMR lines which grows as temperature is 
decreased, in quantitative agreement with recent experi- 
ments by Yamani et ah ^^'^" The Friedel-type spin density 
wave induced by a uncorrelated metallic host is, by con- 
trast, too small and temperature independent and there- 
fore cannot account for the broadening of NMR lines. We 
further proposed an ad hoc model which may account for 
the high temperature satellite peaks in the NMR spec- 
trum observed in Refs. [29lf35. by assuming a set of a 
few percent Cu(lF) 3-1- ions in the chains. The essential 
point is that these exceptional Cu's should produce zero 
polarization; this then accounts quite well for the weight 
and temperature independence of the observed satellite 
peaks. 

Our model assumes pairing interactions in the plane 
alone, and superconductivity is induced in the chains only 
by "proximity coupling" , i.e. a chain-plane electron hop- 
ping. This coupling is then found to cause an imprint 
of the defect-induced magnetic chain correlations on the 
plane itself, consisting of an incommensurate spin den- 
sity wave of quasi-lD character. This signal is an order 
of magnitude smaller than the chain polarization for re- 
alistic parameters but should still be observable. Indeed, 
within our model the broadening of planar NMR lines is 
proportional to that of the chain NMR line, consistent 
with the observation by Yamani et These defect- 

induced ID correlations should enhance the effect of the 
chain bands on the magnetic response of the YBCO sys- 
tem at low energies, particularly for low dopings near the 
Ortho-II and III oxygen concentrations where an integer 
number of chains are filled. Further research into various 
aspects of this response is in progress. 
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